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Abstract 

°: n Q 

As the glass (in molecular fluids [1]) or the jamming (in colloids and grains [2]) transi- 



tions are approached, the dynamics slow down dramatically with no marked structural 



changes. Dynamical heterogeneity (DH) plays a crucial role: structural relaxation oc- 



, y, y. on 



dynamical arrest. Here, we measure in colloids the maximum of a "dynamical suscep- 



curs through correlated rearrangements of particle "blobs" of size 
approaching these transitions, £ grows in glass- formers [a, la], colloids [3, 7], and driven 
granular materials [8] alike, strengthening the analogies between the glass and the jam- 
ming transitions. However, little is known yet on the behavior of DH very close to 

O 

tibility", x* ■> whose growth is usually associated to that of £[9]. x* initially increases 

00 ! 

with volume fraction ip, as in[8], but strikingly drops dramatically very close to jam- 
ming. We show that this unexpected behavior results from the competition between 

x. 

the growth of £ and the reduced particle displacements associated with rearrangements 
in very dense suspensions, unveiling a richer-than-expected scenario. 

PACS numbers: 61.43.Fs Glasses, 64.70.Pf Glass transitions, 82.70.Dd Colloids, 83.80.Hj Suspensions, dis- 
persions, pastes, slurries, colloids 
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The recent observation of a critical-like growth of temporal and spatial dynamical flue 



tnat 10 n S in a 2D athenaal sys te m approaching jamming simile, to that hypothesized 

for glass formers 1 101]. has raised hope that the glass and the jamming transition may be 
unified, calling at the same time for further, tighter experimental verifications. Here, we 
investigate temporal DH in a 3D thermal system, concentrated colloidal suspensions close 
to the maximum packing fraction. Temporal and spatial DH are usually closely related: 
the former can be quantified by a "four-point dynamical susceptibility" %4 (the variance 
of a time-resolved correlation function describing the system relaxation), whose amplitude 



is proportional to £' 



■ay 



12j . Surprisingly, we find that very close to jamming temporal 
and spatial DH decouple: while £ continuously grows with volume fraction, the amplitude 
of temporal fluctuations drops sharply close to the maximum packing fraction. These find- 
ings challenge current scenarios where the slowing down of the dynamics on approaching 
jamming is accompanied by enhanced dynamical fluctuations. 



We study concentrated suspensions of poly vinil chloride (PVC) xenospheres[13l| suspended 
in dioctilphtalate (DOP). The particles are highly polydisperse, with typical diameter « 10 
jum; they behave as slightly deformable hard spheres. The samples are loaded in cells of 
thickness L = 2 mm, vigorously stirred and outgassed to remove air bubbles. The dynamics 
are probed by dynamic light scattering in the highly multiple scattering limit (Diffusing 



Wave Spectroscopy, DWS[l4(), adopting the transmission geometry, with L/£* « 10, where 



is the photon transport mean path. This technique allows the dynamics to be probed on 



length scales as small as a few nm 141]. which match well the restrained motion in tightly 



packed suspensions. A charge-coupled-device (CCD) detector is used to record the speckle 
pattern scattered by the sample. The evolution of the speckle images is quantified by the 
two-time degree of correlation [15( a(q, t, r) = (I p (t)I p (t + r)) p / ((I p (t)) p (I p (t + r))) - 1, 
where I p (t) is the scattered intensity at pixel p and time t and (• • •) is an average over the 
CCD pixels. 

In order to follow the evolution of the dynamics, we calculate g2(t,r) — 1, the two-time 
intensity correlation function obtained by averaging ci(t, r) over a few CCD frames. Due to 
the limited acquisition rate of the CCD, the initial decay of g 2 (t,r) — 1 is not captured; at 
longer delays, a plateau followed by a final relaxation is observed (see Fig. 2a for an example 
of a time-averaged g 2 ), indicative of very slow rearrangements. Figure la shows a typical 
example of the time dependence of tq, the characteristic time of the final relaxation obtained 



by fitting g 2 — 1 to a stretched exponential a(t) exp < — ]t/to(£)] >. Initially, r grows 
nearly linearly with t, as observed in many glassy systems. However, for t > 28000 sec a 
stationary regime is observed, where r exhibits surprisingly large fluctuations but no overall 
increasing trend. A similar behavior is observed for all volume fractions; all data presented 

n 

in the following refer to the stationary regime [16j. We first investigate the (^-dependence of 
the average dynamics, as quantified by To and 3, where 7TT denotes a time average. As shown 
in Fig. lb, Tq continuously increases with tp. Data taken for freshly prepared samples (solid 
circles) can be fitted by a critical law, To ~ 1/l^/Vma* — 11^, with y3 max = 0.752, consistent 
with expectations for highly polydisperse samples 171]. and x = 1.01 ± 0.04, similarly to 



Ref.[8|. We also study samples that have been aged for several days and whose dynamics is 
re-initialized by vigorously stirring and outgassing them (semi-open circles). Their effective 
volume fraction is higher than the nominal one, due to the slight swelling of PVC particles 
suspended in DOP for very long times[13J], allowing to achieve an even tighter packing. 
To compare the dynamics of both fresh and aged samples, we assign an effective volume 
fraction, ip e Q, to the latter so that their average relaxation time falls on the critical-like curve 
determined for the fresh samples. Figure lc shows /?(</?): at the lowest volume fraction, the 
shape of g 2 is slightly stretched (3 = 0.86 < 1), similarly to what observed for correlation 
functions in many glassy systems[l|]. Surprisingly, as (p increases 3 grows above one, finally 
satura ting around 1.3. A similar "compressed" exponential relaxation has been observed in 



single [18|, [19|, |20j and multiple [2 ll. |22J scattering experiments on systems close to jamming, 
usually associated with ultra-slow ballistic motion. 

We quantify the temporal fluctuations of the dynamics by calculating x( T , V 9 ); the relative 
variance of cj, defined by 



X( r ) = X{r, tp) = (ci{t, t) - c 7 (t, r)j /a 2 



(1) 



where the normalization is introduced to account for the (^-dependence of the amplitude 
of the final relaxation of g 2 — 1; data are furthermore corrected for experimental noise 15]. 



The variance introduced above corresponds to the dynamical susceptibility y 4 much studied 



in simulations of glass formers [9|, |ll|, 



Figure 2a shows both the average correlation 



function g 2 — 1 (open circles) and \ (solid circles) for ip = 0.738. The dynamical susceptibility 



exhibits a mar 

glassy systems 



ied pea 

,y,y, 



£ around To, a direct manifestation of DH also found in many other 
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12 



24 1 



Figure 2b shows the height of the peak of the dynamical susceptibility, x*i as a function of 
(p. At the lowest volume fractions, x* increases with (p; the data can be fitted by a critical law 
X* ~ l/l'/'/Vmax - l| y with y = 1.5 ±0.2 (line in Fig. 2b), close to y = 1.70 recently reported 
for driven grains 1 8]. This growing trend is also analogous to that observed in simulations 
of glass formers [5J and colloids [3|, |7| (albeit at lower ip) and has been interpreted as due 
to a growing dynamical length scale on approaching dynamical arrest. At higher volume 
fractions, however, an opposite trend is observed: the amplitude of dynamical fluctuations 
dramatically decreases close to <p max . This striking behavior represents our central result, 
which challenges current views of DH close to dynamical arrest. The unexpected drop of 
dynamical fluctuations very close to jamming is confirmed by the non-mo notonic behavior 
of the width of the temporal distributions of r and (3, shown by the vertical bars in Fig. 
lb-c. The dispersion of both parameters initially increases with (p, but is eventually reduced 
close to <p max , further demonstrating reduced DH. 

We propose that the non-monotonic behavior of x* results from a competition between 



the growth of £ on approaching (p max and the reduced particle disp 



rearrangement events close to jamming, due to tighter packing 
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acement associated with 
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27 



281 ] . Indeed, as 



£ increases, fewer statistically independent dynamical regions are contained in the sample, 



leading to enhanced fluctuations |24(. Conversely, as particle displacement decreases x* is 
reduced, since more events are required to significantly decorrelate the scattered light and 
fluctuations on a time scale ~ tq tend to be averaged out. These competing mechanisms 
should be quite general and should be observable in a variety of systems, provided that DH 
are probed close enough to dynamical arrest. 

We have incorporated these ideas in a simple model for DWS for a dynamically het- 
erogeneous process, significantly extending previous work on the intermittent dynamics of 
foams [291] and gels[30j. The dynamics is assumed to be due to discrete rearrangement events 



of volume £ 3 that occur randomly in space and time; however, in contrast to Ref. 29] we 
assume that several events will be in general necessary to fully decorrelate the phase of 
scattered photons, since in concentrated suspensions the particle displacement associated 
with one single event may be much smaller than the wavelength of the light. The two-time 
field correlation function for a photon crossing the cell along a path of length s may then 
be written as 

sj s) (;,r)=exp[-n s (;,T)Vj]. (2) 



Here, n s (t, r) is the number of events along the path between time t and t + r and o\ is the 
variance of the change of phase of a photon due to one single event, related to the particle 
mean squared displacement associated with such event, a 2 , by cr 2 , ~ 20cr 2 //im 2 [3l|. For a 
totally uncorrelated change of photon phase due to distinct events one has p — 1, while in 
the opposite limit of a perfectly correlated change of phase p = 2[14]. 

We implement our model in Monte Carlo simulations where photon paths are random 
walks on a square lattice with lattice parameter £*. The lattice sites are affected by random 

rearrangement events of volume £ 3 , occurring at a constant rate per unit volume. The 

r Ml 2 

simulated degree of correlation is calculated from c/ )S i m (£, r) = iV" 1 ^2 9i ($■> T ) > where 

the sum is over N s = 200000 photon paths and g[ s is calculated according to (J2J). The two- 
time intensity correlation function, g2, S im(t, r), and its fluctuations, Xsim, are then calculated 
from Cf.sim as for the experiments. We vary the control parameters in the simulation, £ 3 , 
p, and cr 2 ,!^, to reproduce the experimental (^-dependence of (3 and x* ■ Since the particle 
displacement resulting from one rearrangement — and thus <j\ — is expected to decrease as if 



grows, due to tighter particle packing [25 



26 



27 



28(, we choose 1/cr? as the control parameter 



against which simulation results are presented, corresponding to increasing volume fractions 
in Fig. lb-c). Figure 3a shows (3 vs 1/cr 2 ,. The data are obtained using the values of £ 3 shown 
in Fig. 3c; however, we find that /3 depends only very weakly on £ 3 . For large particle 
displacements (small 1/cr?), (3 < 1 in fair agreement with the experimental value at the 
lowest ip. As particle displacements become increasingly restrained, (3 grows and saturates 
at (3 « 1.3, close to the experimental values at the highest if. The saturation value depends 
on the choice of p 33j: here, p = 1.65, showing that the change of phase of a photon due to 
distinct rearrangements is partially correlated. It is unlikely that such a correlation exists 
for events occurring in non-overlapping regions; by contrast, successive events in the same 
region will lead to partially correlated changes of phase when the direction of displacement 
persists during several events. Thus, p = 1.65 indicates intermittent supradiffusive motion, a 

nrv 

behavior close to the ballistic motion reported for many jammed systems [la, H9| 
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22] 



Figure 3b shows the 1/cr 2 dependence of \* ■ The simulations reproduce well both the non- 
monotonic trend and the absolute values of the experimental x* ■ They reproduce also the 
non-monotonic (^-dependence of the dispersion of (3 (bars in Fig. 3a), once again matching 
closely the experimental data. In spite of the drop of dynamical fluctuations close to y? max , 
£ 3 grows steadily with if (Fig. 3c), until rearrangement events span the whole sample, 



corresponding to £ ~ 2000 particle diameters. System-spanning "earthquakes" have been 



reported in simulations of both thermal 34j and athermal 28[ systems close to dynamical 
arrest, but have never been observed experimentally. 

Our results show that the behavior of dynamical heterogeneity very close to the jamming 
transition is much richer and complex than expected. On the one hand, the growth of 
dynamically correlated regions is limited only by the system size, implying that confinement 



effects, usually observed on the scale of tens of particles at most |35| , should become relevant 
macroscopically. On the other hand, temporal fluctuations of the dynamics are suppressed, 
challenging current views of DH close to jamming and calling for new theories. 
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FIG. 1: Time and volume fraction dependence of the relaxation of concentrated suspen- 
sions, a) Time dependence of the final relaxation time obtained by fitting the two-time intensity 
autocorrelation function by a stretched exponential, giit, T ) = o,(t) exp[— (t/to(£))^*'], for a sample 
at (p = 0.664. After an initial aging regime where the dynamics slows down, the system reaches a 
dynamically heterogeneous stationary state, where tq fluctuates significantly without any overall 
growing trend (the time origin is taken at the end of the sample loading and outgassing). The 
lower panels show the volume fraction dependence of the relaxation time (b) and the stretching 
exponent (c) averaged over time in the stationary regime. The solid symbols refer to freshly pre- 
pared samples, the semi-filled circles to old samples whose dynamics has been re-initialized. The 
bars are the standard deviations of the distributions over time of tq and (3 in the stationary regime. 
The solid line in b) is a critical law fit to the growth of tq for the fresh samples, yielding a critical 
exponent x = 1.01 ± 0.04 and a critical packing fraction <£> max = 0.752, indicated by the dashed 
line here and in c). The solid line in c) is a guide to the eyes. 



10 



0.3 



I I llllllj 1 I llllllj 1 I Mill 1 I I llllll TTT 

o Oct) 



0.02 




N 



0.01 



0.00 



0.00 



0.60 0.65 0.70 0.75 



<p,<p 



err 



Figure 2; Manuscript NPHYS-2007-05-00481D 
P. Ballesta, A. Duri, and L. Cipelletti 



FIG. 2: Dynamical susceptibility, a) Average correlation function 52 ( r ) — 1 (open symbols and 
left axis) and dynamical susceptibility (solid symbols and right axis), for a sample at ip = 0.738. 
b) Volume fraction dependence of the height of the peak of the dynamical susceptibility (same 
symbols as in Figs. lb-c). The solid line is a critical-law fit to the initial growth of %*, yielding 
an exponent y = 1.5 ± 0.2. Note the unexpected drop of x* near the maximum packing fraction, 
shown by the dashed line. 
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FIG. 3: Simulations. Time- averaged stretching exponent (5 (a) and peak of the dynamical sus- 
ceptibility x* (b) as a function of 1/c?, as obtained from the simulations described in the text (note 
that l/<7? increases with ip). In a), the bars indicate the standard deviation of the distribution of 
stretching exponents. The volume £ 3 = (£/^*) 3 of the rearranged regions used as an input in the 
simulations in order to reproduce the cp dependence of (5 and x* is shown in (c). £ saturates at the 
value of the smallest dimension of the scattering cell. 
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